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ABSTRACT 


An analysis has been presented for determining the 
dynamic response of an articulated non-linear landing gear 
at touch-down. The analysis with the aircraft modeled as a 
rigid lumped mass in heave has been extended to a heave- 
pitch model with an articulated nose landing gear and 
telescopic main landing gears. In developing the system 
equations oil compressibility effects, kinematics of the 
articulation of the gear, drag loads on the tyre and 
pneumatic, hydraulic and frictional forces in the shock 
strut have been included. A parametric study of the orifice 
discharge coefficient, polytropic index of air compression 
process, orifice diameter, initial air volume and initial 
air pressure on the heave model and the pitching moment of 
inertia and the landing gear layout configuration on the 
heave- pitch model of the aircraft have been carried out. 
Numerical results have been obtained for a light fighter 
trainer aircraft by numerically integrating the coupled non- 
linear differential equations of motion. 
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CHAPTER 1 


INTRODUCTION 

In order to cater to the requirements of landing impact 
and ground roll, the landing gear of an aircraft 
functionally consists of a wheel unit and a shock absorber 
unit. In the past, the shock absorbing characteristics of an 
aircraft were developed by means of extensive trial and error 
drop testing. However, the drop testing does not properly 
simulate the aerodynamic forces on the aircraft. Moreover, the 
expense and time required for such an approach emphasises the 
need for suitable theoretical methods and appropriate 
computer simulation for the analysis of the landing gear 
behavior. Such a rational method finds application in : 

-determining the behavior of a given landing gear to 
varying impact conditions. 

-developing new landing gear configurations to suit 
given impact condition. 

-determination of wheel spin-up and spring-up loads which 
takes into account the shock absorbing characteristics. 

-determination of dynamic loads in flexible aircraft 
structures during landing. 

-estimation of airframe structural fatigue and 
fatigue of undercarriage load attachment points. 

Apart from the effect of vibration on airframes, it is 



also necessary to consider the influence of vibration 
environment on crew and passengers. In some cases, the cockpit 

vibration in the vertical plane can degrade the pilot’s ability 

to read the instruments accurately and perform correct control 

operations. Simulation studies suggest that simultaneous 

exposure to vertical acceleration greater than 0 . 6g at 3 Hz and 

horizontal acceleration greater than O.lg at 2 Hz may cause 

some loss of visual acuity [1]. It is known that cockpit 

vibration levels are greater than that experienced at aircraft 

mass center due to the pitching motion. Thus pilots will be 

more influenced by vibration than passengers. With regards 

to the effect of critical frequencies on human beings, 

human response to both harmonic and random vibrations is 

most critical between 1 Hz and 5 Hz and it is in this range 

that most air-frame vibrations are excited [2]. 

Since many aspects of landing impact problems are so 
intimately connected with the mechanics of landing gear, the 
subject of landing gear behavior has received analytical 
treatment at various times. Many of the earlier investigations 
in order to reduce the mathematical complexity of the 
analysis, were limited to consideration of highly simplified 
linear systems, which have little relation to the actual 
landing gears. Some of the later works, consider with varying 
degrees of simplification, more realistic non-linear systems. 

Taxiing induced vibrations in large aircraft due to 
run-way unevenness has also been recognized as a 

significant factor in metal fatigue damage, dynamic stressing 
of airframes as well as crew and passenger discomfort. Such 
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taxiing loads applied to an airframe should not be ignored 
since their cumulative effects represents a considerable 
proportion of fatigue damage. The fatigue life of the 
undercarriage at its point of attachment with airframe 
as well as associated structural components should also be 
considered as an Important design case. Many investigators 
have studied the effects of runway induced vibrations on 
landing gear. 

Based on geometrical disposition, landing gears are 
classified as telescopic or non-articulat ed , levered or 
articulated and semi-levered or semi-articulated [3]. Though 
there are a number of aircrafts with articulated and semi- 
articulated landing gears, most of the published literature 
deals only with the behaviour of telescopic landing gears. the 
present work analyses heave and pitch models of an aircraft 
with articulated landing gears. 

l.i REVIEW OF PAST WORK 

Since the dynamic conditions present during landing 
impact are often critical, from very early days considerable 
emphasis was given to the accurate prediction of landing gear 
behavior in the landing run. 

Schlaefke [4] analyzed a telescopic landing gear 
in which the spring was considered to be linear and for the 
first time assumed the damping to be proportional to the 
square of the spring compression velocity. From the study, he 
drew a disputable conclusion that the realization of quadratic 
damping law would offer no advantage over linear damping. 



Milwitzky and Cook [5] made a theoretical study o£ 
behavior of a telescopic oleo-pneumatic landing gear during 
the process of landing impact. The analysis treated the motions 
of landing gear prior to and subsequent to the beginning 
of shock strut deflection. In the first phase of landing 
impact, a single degree of freedom system was solved in 
order to determine the conditions at the instant of initial 
shock strut deflection, after which a two degrees of freedom 
system was solved. The analysis considered such factors 
as hydraulic (velocity squared) resistance of orifice, the 
forces due to air compression and internal friction in the 
shock strut, non- linear force-deflection characterist ics of 
the tyre, wing lift, the inclination of landing gear and the 
effect of wheel spin-up and drag loas. In addition to this 
exact treatment, investigations were also made to determine the 
extent to which the basic equations of motion could be be 
simplified and yet yield acceptable results. Most of the 
specific solutions were obtained by "linear procedure” which 
assumes the changes in the variables to be linear over finite 
time intervals. A few of the specific solutions were also 
obtained by "quadratic procedure" and the generalised 
solutions of equations of motions were obtained with the 
well known Runge- Kutta method. 

An analog simulation of landing impact and initial 
ground run was performed by Uahi [6]. He included in his work 
the internal friction in shock strut, wing lift, inclination of 
landing gear and effects of wheel spin-up and drag loads. The 
discharge coefficient was treated as function of Reynold’s 
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number, orifice shape and orifice orientation. An iterative 
procedure was adopted for the solution. His later paper [7] 
included the effects of hydraulic fluid compressibility, gas 
solubility, and its entrainment in oil, and the varying nature 
of polytropic exponent during compression stroke. The model 
chosen by Uahi was a heave model in which the pitching of the 
aircraft was not accounted for. 

Nightingale [8] described an approximate me thod for 
predicting landing gear performance from a 'step-by-step* 
calculation by assuming a load-stroke curve for an articulated 
landing gear. Equations were developed £ 0r a telescopic landing 
gear with a vertical strut and later modified to include 
telescopic landing gears with inclined strut and levered 
landing gear. A simple model was considered for the analysis. 
The aircraft was modeled as a single mass, supported by an 
oleo-pneumatic shock strut and a tyre in series. The tyre was 
taken to be an undamped linear spring and the landing was 
assumed to be fully airborne. The effects of surface friction 
and the mass of wheel were not considered in the analysis. 

Others also have analysed the behavior of articulated 
landing gear as extension of telescopic landing gear. In all 
these works, a single degree of freedom system was 
solved neglecting the unsprung mass. Mechanical advantage and 
velocity ratio were obtained from layout drawings of the gear 
and were used for obtaining landing impact transients. The 
kinematic relation between wheel axle travel and shock strut 
deflection were not introduced in these mathematical models. 

Jayarami Reddy et al [9] developed dynamic equations 
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with kinematic relations as a part of analysis for an 
articulated landing geai . He consideted in his work velocity 

squared damping, seal and journal friction and treated the tyre 
as an undamped non-linear spring. On impact only the tyre was 
assumed to compress and hence a single degree of freedom system 
was solved. Once the force in the shock-strut exceeded its pre- 
load, the shock-strut became active and a two degrees of 
freedom system was solved. A fourth order Runge-Kutta method 
was used for solving the coupled non-linear equations of 
motion. A semi- levered suspension type landing gear was also 
analysed by the same authors [10]. This investigation included 
the effects of oil compressibility and considered the wheel 
spin-up as a function of slip ratio. Both these works assumed 
only a heave model in which the pitching motion of the aircraft 
was not considered. 

1.2 SCOPE OF THE PRESENT WORK 

The present work develops an analysis for the 
determination of the response of an aircraft at touch-down. 
Initially the aircraft is idealised to a single rigid mass 
supported on a levered landing gear. The hydraulic, pneumatic 
and friction forces of the shock strut and the linkage 
dynamics of the gear have been taken into account for the 
study. Other factors included are the wing lift forces and the 
drag loads on the wheel during its spin-up. 

Due to the complicated nature of the thermodynamic 
process occurring in the shock strut at impact no consensus 
exist on the exact value of the polytropic index of air 
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compression at these condition. The coefficient of orifice 
discharge and the area of orifice also varies with the shock 
strut stroke due to the presence of metering pin and is 
diffetent for compression stroke and expansion stroke. 
Therefore the variation in the dynamic response of aircraft to 
different values of these parameters are investigated. 

To make this model more realistic the pitching motion 
of the aircraft is Included in the second model. The aircraft 
considered is assumed to possess a levered nose gear and 
telescopic main landing gear. Since the main gears are 
usually located fairly close to the nodal points of the 
fundamental bending mode of the wing, the landing gear 
performance is not significantly affected by the elastic 
deformtion of the wing. In the case of large supersonic 
transports which have large slender bodies, elastic modes in 
the longitudinal direction may also be considerable For 
smaller aircrafts the airframe may be assumed to be rigid 
without significant loss of accuracy. The flexible modes may 
also be incorporated in this study by the simultaneous 
solution of the equations of motion of the landing gear 
coupled with those equations representing the additional 
degrees of freedom of the structure. 

Rolling and yawing of the aircraft is not included in 
developing the system response. The pilot’s input to the 
system is also not considered because of lack of reliable 
information regarding pilot’s action at touch down. 

The effects of the variations in the pitching radius of 
gyration and the location of landing gear with reference to 
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the center of gravity of the aircraft is also studied. 

By changing the parameters described above a spectrum 
of responses is developed. In all these cases the following 
quantities were determined : 

a. translation, velocity and acceleration of the wheels. 

b. vertical translation, velocity and acceleration of the 
center of gravity of the aircraft. 

c. pitching of the aircraft about its center of gravity. 
The dynamic response so obtained can be used for 

determining the load factors, ground clearance and other 
standards specified by the international agencies to determine 
whether any damage has occurred. 

1.3 LAYOUT OF THE THESIS 

The general layout of the thesis is presented below. 
Idealization of the aircraft and the details of the 
development of system equations are dealt in chapter 2. A 
parametric study on the landing gears of the aircraft and its 
effects on the aircraft response during landing impact are 
covered in chapter 3. Chapter 4 presents the conclusions drawn 
from the study and suggestion for further investigations in 
this area. 
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CHAPTER - 2 


MODELING OF THE LANDING GEAR 

The accuracy with which the behaviour o£ any dynamic 
system is predicted depends to a large extent on the model 
chosen to idealize the physical system. The present work 
analyses the response of an aircraft with articulated nose gear 
and telescopic main gears. Two models have been considered for 
the analysis. First, a heave model is studied for a levered 
gear with linkage dynamics included as a part of the analysis. 
This is extended to a heave-pitch model for the aircraft which 
Includes the full kinematics of the articulated nose gear and 
conventional telescopic main gears. The landing gears are 
assumed to be infinitely rigid in bending and the external lift 
forces are assumed to act throughout the impact and 
ground-roll. The analysis also includes the oil compressibility 
effects, pneumatic, hydraulic and friction forces in the shock 
strut . 

2.1 MECHANICS OF THE LEVERED LANDING GEAR 

The type of articulated landing gear analyzed is shown 
in fig 2.1. The shock strut (ED) of the gear is hinged at both 
ends and hence experiences only axial loads. The trailing arm 
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FIG. 21 ARTICULATED LANDING GEAR 
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(AC) has the wheel axle at its lower end (C). The housing of 
the shock strut is mounted on the aircraft which has a fork 
hinge (A) at its bottom for receiving the upperend of the 
trailing arm and fork end (E) at its top for receiving the 
cylinder eye end of the shock strut. The piston eye end (D) is 
hinged to the trailing arm through a rigid link (BD). The 
hinges at A and E are thus fixed with references to the 
airframe. 

Figure 2.2a shows the geometric layout of the landing 
gear in the fully extended condition when the aircraft is 
airborne. During the process of landing, the trailing arm moves 
from its original position AC to the deflected position AC’ in 
an arc of circle with A as center, shifting the axis of the 
shock strut from ED to ED' (fig 2.2b). The relative motion 
between the aircraft mass and the wheel compresses the shock 
strut (ED) producing the necessary resistance in the form of 
pneumatic, hydraulic and friction forces. 

The inclination of the trailing arm with the vertical <9 t 
at any instant of time t is 


6 

t 



. . . ( 2 . 1 ) 


where 

A 4 ..A tf = geometrical parameters describing the configuration 
of the landing gear in the fully extended condition 
z * vertical displacement of the aircraft mass 
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vertical travel of wheel axle 


z 

a 


The extended length of the shock strut s^ 
from the geometry. Thus, 


is obtained 




r * * * 1 

S o = | [k v Since o +<5)-Aj + [k 9 CosCd o +<5)+A B )] + [y B -3^ 3 j 


. . .( 2 . 2 ) 


where 


y . . .y = 


© 


4X2 


K-O 

-i 

Tan (A 4 / A a ] 

lateral distance of the shock strut joints A...E 
from the center line of the aircraft 

inclination of the trailing arm from the vertical in 
the airborne configuration of the landing gear 


The compressed length of the shock strut s^ at any 
instant of time t can be obtained from equation (2.2) by using 
the current value 6^ for the inclination of the trailing arm in 
place of & o 


4/1 

f * 2 * 1 

s t = 4 [A^ Sin(0 t +<5)-A 4 J + [A 7 Cos(O t -h5)+A B )3 + j 

...(2.3) 
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The shock strut stroke is therefore s = s . - s. ...(2.4) 

Differentiating equation (2.4), the telescoping velocity 
of the shock strut is obtained. Thus, 


A T (z t - z 2 )j A e Sin(€? t +£) + A i Cos(© t +<5) j- 


A Sin & s 

a t t 


. .(2.5) 


From the geometry, the horizontal travel x of the wheel 


axle is 


x = A ( Sin & - Sin & ) 
a t o 


. . . ( 2 . 6 ) 


Differentiation of the above expression gives the 
horizontal velocity and acceleration of the wheel as 

x = q ± ( z a - z 2 ) ...(2.7) 


x = q 4 (z 4 - z a ) - <l 2 (2 i - z a ) 


where 




Cot & 


l/A a Sin a 


( 2 . 8 ) 


(2.9) 


a. a ANALYSIS OF THE HEAVE MODEL 

The dynamic system considered for the analysis is shown 
in figure 2 . 3 .Oleo-pneumatic struts are inflated to some finite 
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air pressure in the fully extended position. The strut doesn't 
begin to deflect in an impact until sufficient force Is 
developed to overcome the preload F a& due to the extended air 
pressure and internal friction. Since the struts are 
comparatively rigid in bending and compression , initially the 
shock strut is locked and the tyre alone compresses with the 
system behaving as a single degree of freedom system. Once, the 
preload is exceeded, the shock strut becomes operational and 
the system is converted to one with two degrees of freedom. 

2.2.1 Equations of Motion with Shock Strut Locked 

In the locked condition of the shock strut, the aircraft 
and the wheel moves as a whole and hence the system is modeled 
as a single degree of freedom system. 

Equation of motion for this model in the vertical 
direction is 


U . . 

-2+L-U+F =0 .. .(2.10) 

£ 

where 

U = portion of aircraft weight acting on the landing 

gear 

L = aerodynamic lift force on the aircraft 

F » vertical ground reaction on the tyre 

v«gr 

2 = vertical translation of the aircraft and the wheel 

mass 

From the force equilibrium of the sprung mass in the 
vertical direction, the vertical component of the force in the 
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trailing arm F is obtained as 

VOL 

F = U - L 

VO 1 


u . 

4 


z 


fi 


. . .( 2 . 11 ) 


where is the weight of the upper sprung mass 

Zn the locked configuration of the landing gear the 
horizontal component of the force in the trailing arm F^ can 
be determined from the force equilibrium of the wheel as 


F = - F 

Kol hg 


. . .( 2 . 12 ) 


where F L is the horizontal ground reaction on the wheel due to 

Kg 

ground friction. 

The inclination of the axis of the shock strut 4 > with 


the vertical at any instant of time t is 




-4 

Tan 


A Sin(0 +6) - A 

? t 4 


A Cos(<9 +6) + A 

*? 1 B 


. . .(2.13) 


The force F in the shock strut is obtained from the 

a 

balance of moments about the hinge A of the trailing arm. Thus, 


where 


F 

* 


F - 

VOL 


* Ha 


q 


* 


q 


B 


quantity defined in equation (2.9) 


% 


A ? [ l a Sin(e t +6) - l K Coa(S t +<5) ] 

A Sin & 

a t 


. . . (2.14) 
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1 = Cos $ 

3 *■ 

1 = Sin A 

K 

The force F in the shock strut is monitored 

a 

continuously. The system transits from one to two degrees of 
freedom system when the force in the shock strut exceeds its 
preload F . The state at which F = F gives the initial 
condit ion for the solution of the two degrees of freedom 
system . 

2.2.2 Equations of Notion With shock Strut Unlocked 

Once the shock strut becomes operational , the sprung 
mass and the unsprung mass displace separately therefore the 
system is modeled as one with two degrees of freedom. 


Vertical equilibrium of the sprung mass is 

U 

* .. ...(2.16) 
F + L - U + — z = 0 

V«& 1 g 4 

Equation of motion of the unsprung mass in the vertical 
and horizontal directions are 



IS 


. . .( 2 . 19 ) 



where 


q„ 


U 

a 

T q * 


¥ . a 

= “ l ^ V + ^ 3 


( 2 . 20 ) 

= - r— afz-z'l+F 1 

4 

q^ 4 q a quantities defined in equation (2.9) 

From the moment equilibrium of the trailing arm about 
the hinge A, the vertical component of the force in the 


trailing arm F can be obtained as 

v<a. 

F = q F + q F 

VOL B • 4 ha 

Eliminating F from equations (2.19) and (2.21) 

FiOL 

on simplification 

F * B + B (**-«.) 

vol a a * 1' 


.(2.21) 

yields 

.(2.23) 


where 


B 

s 

q F 

+ q.q. 

ft 


fir ft 

i 4 

B 

a 

q q^ 


a 


A Tl 


Let 




B 

s 

B + 

L - U 

* 


& 

4 

B 

35 

B ♦ 

u - 

6 


a 

a 

B 

a 

B + 

W / g 

<s 


a 

4 

B 

55 

B + 

U /g 

9 


» 

a 

B 

ft 

= 

B* - 

ft 

B B 
a ? 




(2.24) 


Substitution of equations (2.22) and (2.24) in the 
equations of motion (2.16) and (2.17) yields on simplification 

. . .(2.25) 


z M BB - B B 3/B 

A * 9*4 ft fir f ft 
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Z = [ B B - B B ]/B 

The coupled non-linear differential equations 
(2.26) are the system equations of motion which 
numerically by standard routines. 


. . . (2.26) 
(2.25) and 
are solved 


2.3 ANALYSIS OF THE HEAVE-PITCH MODEL 

The heave model considered in section 2.2 is improved to 
a heave pitch model by considering the operations of the nose 
and the main wheels separately so as to account for the 
pitching motion of the aircraft. Figure 2.4 shows the model 
used for the analysis. 

The combination of aircraft and landing gear consists of 
four degrees of freedom due to the rigid body translation and 
rotation of the upper sprung mass and the vertical displacement 
of the lower unsprung masses of nose and main gears. All motion 
is assumed to be restricted to the vertical plane through the 
center of gravity of the aircraft . 

According to the force in the shock strut in the main 
and nose gear, four distinct possibilities exists. 

1. Both main and nose gears are locked. 

2. Nose gear is locked and the main gear is unlocked. 

3. Main gear is locked and the nose gear is unlocked. 

4. Both main and nose gears are unlocked. 

These four cases are considered one by one in the 
following sections for developing the equations of motion. 
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FIG. 2-4 THE HEAVE- PITCH MODEL 
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2.3.1 Equations of notion with Main Gear and Nose Gear Locked 

Uhen the forces in the shock strut are less than their 
preloads, only the two tyres compress and hence the system has 
two degrees of freedom with the unsprung mass moving along with 
the airframe. 

The equation of motion of the system in the vertical 
direction is 


where 


U . . 

-z + L- U+ F + F =0 

g 4 vq vof 

to r» 


. . . (2.27) 


U 

U 


U 


U, 


z 

A 


= + W ■* U, (total weight of the aircraft) 

= weight of the sprung mass 
= unsprung weight of the nose wheel 

- unsprung weight of the main wheel 
= lift acting on the aircraft 

- vertical displacement of the sprung mass 


F 


v «. 


F 




TO 



vertical ground reaction on the wheels 


subscript m for the main gear 
subscrpit n for the nose gear 

Since the shock struts are locked the displacements z 

ft 

of the noae wheel and z of the main wheel are obtained from 
the geometry. Thus, 


z 5 

ft 


z - 1 z 

4 n ft 


Z « 

4 


£ * 1 2 
A m & 


. . . (2.28) 
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where 


distance between aircraft c.g and the nose wheel 
distance between aircraft c.g and the main wheel 
angle of pitching of the aircraft 


Differentiating twice, the acceleration of the nose 
wheel and the main wheel are obtained as 


2 * Z - 1 Z 

a a r. a 

z * z 1 z 
* a m a 


( 2 . 29 ) 


Equation of motion of pitching of the aircraft is 


I z - F 1 - L 1 + F 1 + F w (h 

yy ft vg r» oc vj to % 


z ) + F t (h - z ) = 0 

A A ' 

fS 


( 2 . 30 ) 


where 


I = pitching moment of inertia of the aircraft about its 

yy 

c.g. 

1 * distance of the aerodynamic center of the aircraft 

AO 

measured in the forward direction from the aircraft 


c-6- 

F^ = horizontal ground reaction on tyre 

h = height of the aircraft c.g. from the ground level 

The horizontal equilibrium equation of the nose wheel is 


F 


ha 



. . .( 2 . 31 ) 


where is the horizontal component of the force in the 
trailing arm. 

Considering the sprung mass separately the vertical and 
the rotational equlibrlum equations are 
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0 


. . .(2.32) 


U 


- z + L - U + F + F 

g 1 i va a 


I z - F 

yy ft va 


L 1 + F 


1 + F. 

m Kg 




- F CA Cos & 

ha ft O 

n 


1 z ) = 0 

pi a 


. . . (2.33) 


where 

F = force in the shock strut 

A 

A o = dlstancs of the point of attachment of the main e.ar 

from the axle of the main wheel 
Eliminating F from equations (2.32) and (2.33), the 
vertical component of force in the trailing arm is 


F 

va 

n 


u 


A 




A CosO -1 z 
a o n a 


1 


} 


. . .(2.34) 


The moment equilibrium of the forces in the trailing arm 
about the hinge A gives the force Fan in the shock strut of the 
nose gear. Thus 


F - F q 

va ha i 

F^ = — ...(2.35) 

% 

where q^, are terms defined in equations (2.9) and (2.15) 

F 

The force s in the shock strut of the main gear is 

m 

p 

obtained by substituting for va from equation (2.34) in 

TV 

equation (2.32) 
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The system has one increased degree of freedom when the 
force in any of the shock strut exceeds its preload. The state 
at which F exceeds F gives the initial condition for the new 

ft mO 

system. 


2.3.2 Equations of Motion with the Nose Gear Locked and the 
Main Gear Unlocked 

The three degrees of freedom system in this case is 
attained when the main shock strut alone is operating. 


The vertical displacements of the attachment points z 

r» 

and z ffi of the nose gear and the main gear respectively due to 
translation and rotation of the airframe are 


z = z - ( 1 + A Sin 6 )z 

n 4 n a o s 

z = z + 1 Z 

m * tw 2 


. . .(2.36) 


Differentiating once the corresponding velocities are 
z 


= z - ( 1 + A Sin 6 )z 

4 n a oft 


z = z + 1 z 

m 4 m ft 


. . .(2.37) 


The shock strut stroke s and the strut telescoping 


velocity s ^of the main gear are given by 

s = z - z 

m m 4 

S * z -z 

rn m <4 

F 

Since F * F (s,s), the force a 

a a m 

strut can be determined. 


. . .(2.38) 
of the main shock 
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The airframe mass and the mass of the nose wheel acts as 


a single unit. Hence the vertical equilibrium equation of the 
airframe is 


U . . 

— z + L- U+ F + F =0 ...(2.39) 

g 4 vg & 

n m 

where U = U + U 
a i a 


Equation of motion of rotation of the airframe is 
I z - F 1 - L 1 + F 1 +F (A + 1 z ) + F (h -z ) = 0 

yy a v« r. oe • m Kg O m ft Kg 4 

r» tw m r» 

. . .(2.40) 

Acceleration z of the nose wheel in the vertical 

a 

direction is 

z - z - 1 z ...(2.41) 

ft 4 r» a 

Vertical equilibrium equation of the main wheel is 


U 

— z - U + F -F =0 

£4 4 vg o 

w m 


. . . (2.42) 


The equations of motion (2.39) to (2.42) are solved 

simultaneously. The force in the nose gear shock strut is given 

F F 

by equation (2.40), where vo is obtained from (2.32) and h* 

r» r* 

F 

from (2.31). By successively calculating the values of a , the 

1"'» 

F F 

state at which a > a© can be estimated leading to the 

r* 

initial condition for four degrees of freedom system. 
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a. 3. 3 Equations of Mot, ion with the Main Goar Locked and the 
Nose Gear Unlocked 

Uhen the shock strut in the nose gear is operating, the 
linkage geometry relations developed in equations (2.1) to 
(2.7) for the heave model are valid except that z^ and z^ 
correspond to z and z respectively. Thus, 


9 

t 


Cos 


A - (z - z ) 

<s r. a 


. . . (2.43) 


1X2 


S 


f 2 2 2 1 

= -j [A Sin(f9 +d)-A ] + [ A Cos(£ *5 )+A ) ] + [y -JL 3 \ 

[ 7 ° a u *• j 


(2.44) 


i/2 


* t » | [A^ Sin(© t +6)-A A 3 + [A^ Cos^+6 )+A s ) ] + [y^-y^,] J 


s = s - s 

F, o t 


. . . (2.45) 
. . .(2.46) 


s = 


A (z - z )T A Sin(S +«5) + A Cos(0 +<5) 1 

?r, a^= t i ij 


A Sin 9 s 
a t t 


. . (2.47) 


x = A ( Sin 9 - Sin 9 ) . . .(2.48) 

3 to 


x = q 4 (z n - z a ) 


. . . (2.49) 


X = q.(z 

£ D 


z 3 ) ~ 


z a > 
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where the terms used are the the same as those defined 


in equations (2.1) to (2.7) and (2.36) and (2.37). 

• p 

Using s and s , the force a of the nose gear can be 

r* r* r» 

determined. 

The equation of motion of the sprung mass in the 
vertical direction is 


U. 


z + L-U+F + F =0 

i b vcf va 

m n 


(2.50) 


where U = U + W 

b t 4 


Rotational equilibrium equation for the airframe is 
Iz-F 1 -LI + F 1 + E (h - z ) 

yy ft va r» ac vg ro hflf 4 


- F l (A Cos & ~ 1 z_) = 0 

ha ft i r» Z 

n 


. . (2.51) 


Vertical equilibrium equation of the nose wheel is 


U 


— z 
ft * 


U + F - F =0 

a vet va 

n n 


.. .(2.52) 


Since the main wheel and the airframe acts as a single 
mass, the acceleration of the main mass is 


z = z + 1 z 

4 4 m ft 


(2.53) 


The horizontal equilibrium equation of the nose wheel is 
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. . . (2.54) 


U 

3 

F a X - F 

ha 9 h9 

r. r 


Substituting for x from equation (2.49) in equation 
(2.54) yields 


F, = q ( z - z )+q 

ho. a n a * 


where 


U 

a 

r q * 


- [ 


w 

s 


(*-*)+ E ] 

n a K9 


. . .(2.55) 


. . .(2.56) 


Balance of moments of the forces in the trailing arm 

F 

about the hinge A yields the expression for va as in equation 

Yi 

(2.22). Thus, 


F 

n/d 


n 


B + B (z 


z ) 
a 


. . .(2.57) 


where 

B 

z 

B 

a 

F F 

Substituting for from equation (2.57) and ha from 

Ft Pi 

equation (2.55) in the equations of motion (2.50) to (2.53) and 
rearranging in a matrix form yields 

CP] i Z l + i R i - i 0 l ...(2.59) 






(2.58) 
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where the elements of the p matrix are : 


p 

44 


». ♦ «„ / « 

P 

42 

ss 

* \ sin V 

P 

42 

s 

-B 

2 

P 

44 

3S 

0 

P 

24 

« 

B 1 + q (A Cos e - 1 z ) 

2 n 2 2 t n 2 

P 

22 

SB 

V-'yy-'a^a 

P 

22 

S 

-[B aV q a (A fl C0S e t - l n z£ 

P 

24 

= 

0 

P 

24 

SB 

B 2 

P 

22 

= 

-B (1 + A Sin e ) 

a n a o 

P 

22 

SB 

-CB a + U 9 / ft) 

P 

24 

SB 

0 

P 

44 

SB 

1 

P 

42 

SB 

1 

TO 

P 

42 

= 

0 

P 

44 

SB 

-1 


Cos e -1 z )(1 +A Sine 

t nan ni o 


) 


and the elements of the R matrix are : 


R 

s 

B + L - U 

. *• 

+ v<* 

44 


a b 

Tft 

R 

SB 

B 1 + L 1 

F . 

• vg 1 

24 


ft r» ao 

w i 

R 

s 

B - + 

w 

24 


a r, 

ft 

R 

44 

S 

0 



hg ( h- z ) + q (A Cos e 

m 4 4 S t 


1 z ) 

n a 


The force 


F 

a 


m 


in the main shock strut 


can be obtained 
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from the equilibrium equations of the sprung mass in the 
vertical direction. 


F = U 

• 4 


L - F 


W . . 

4 

2 

8 * 


. . . (2.60) 


The system is converted to four degrees of freedom 
system when the preload in the main shock strut is exceeded. 


2.3.4 Equations of Motion with Both Main Gear and Nose Gear 
Unlocked 

Uhen both the shock struts are operating, the forces in 
them can be determined directly. The strokes in the shock 
struts are 


2 - z 


s = s - a 
not 


(2.61) 


The strut telescoping velocities of the main and nose 
gear are obtained from equations (2.38) and (2.47) 


respectively. Since F^ ■ F a (s ; s), the force in the main gear 

F F 

and the nose gear a and a can be determined. 

m r» 

The equations of motion of the sprung and unspung masses 


are 


U 


-iz + L- U+ F + F =0 

£ A A vet A 


(2.62) 


I z ~ F 1 

yy 2 va n 
n 


L 1 + F 1 + F w (A-S*+l z ) 

etc cm her O ?v* m 2 

m m 


- F l (A Cos & - 1 z ) = 0 

ho. a t r, 2 

r> 


(2.63) 
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0 


. . .(2.64) 


U 


g 


- U + F 
a vcf 


- F 


U 


g 


- U + F 


- F 


* 0 


(2.65) 


v s_ 


F F 

Substituting for from equation (2.55) and from 

equation (2.57), the equations of motion (2.62) to (2.65) can 

be written in the form 

t PI *Z* + *R>-<0* ...(2.66) 

where the elements of the p matrix are : 


p 

11 

= 

B a + / g 

p 

18 

* 

-B (1 + A Sin 6 ) 

a n a 6 

P 

IB 

= 


P 

as 

0 

P 

81 

= 

B . l „ + 1 a (A . coa V 

P 

88 

as 

- B . l „<V \ Sin a o 

P 88 

as 

-t B 1 + q_ (A COS & 
*■ a i> a a 

P B4 

= 

0 

P 

Si 

as 

B 

a 

p 

82 

s 

-B a (l n + A fl Sin d o ) 

P 

88 

s 

-(B a + U a / g) 

P 84 

as 

0 

P 41 

= 

0 

P 

4ft 

as 

0 

P 

4ft 

« 

0 

P 

44 

= 

U /g 

4 


yy 
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The el ements of the R matrux are : 


R 

=B 

+ L - 

u 


ii 

a 


4. 


R 

=B 

1 + L 

1 

_ 

BA 

a 

r» 

•CL C 


R 

=B 

F 

- vg 

+ 

u 

a a 

a 

F» 

i 

R 

—I 

w + 

F 


4A 

4 

a 





m 



'vg ] 
m 


2.4 FORCES IN THE SHOCK STRUT 

The total force In the shock strut is due to hydraulic 
resistance, air compression and internal friction. As the strut 
compresses, the hydraulic fluid is forced through the orifice 
to the the lower chamber providing the non-linear damping. Nhen 
the fluid enters the lover chamber, it pushes the separator 
down and compressing the air in the processes. The polytropic 
compression of air provides the non-linear springing. 

2.4.1 Pneumatic force 

The pneumatic force in the shock strut is determined by 
the initial air pressure, initial air and oil volumes, 
pneumatic area oil, bulk modulus and the air compression 
processes. 

The pneumatic force in a shock strut is given by 


33 


where 


P A 

a a 

a 


A a + V 


r^} 


. .(2.67) 


K 

P c 

n 


initial air volume 
initial oil volume 
initial air pressure 
instantaneous air pressure 

polytropic exponent of the air compression processes 
bulk modulus of oil 


2.4.2 Hydraulic force 

The hydraulic resistance F h due to the telescoping of 


the strut is given by 


P A. 


F k= 


s 


2 < c d A . 5 


( 2 . 68 ) 


where 



P 

ss 

mass density of fluid 

\ 

ss 

hydraulic area 

C d 

ss 

coefficient of discharge 

A n 

ss 

net orifice area 

a 

ss 

strut telescoping velocity 
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2. 4. 3 Friction f orcos 


The journal friction F j can be written ae 


. . . ( 2 . f 9) 


where 


u = coefficient of friction of upper bearing (attached 

to the inner cylinder) on the shock strut 


u = coefficient of friction of lower bearing (attached 

to the outer cylinder) on the shock strut 
= normal force on the upper bearing 
= normal force on the lower bearing 
s = strut telescoping velocity 

From the balance of moments 



where 

F 

Ft 


1 + & 
i 




and F = 


1+1 
i 2 


1 + S 

v. l 


so that 


F tJ ■ 


f \ - * 


1 + s 

i 


( P + V. ) + V 




-(2.70) 


normal load on the shock strut 

axial distance between upper and lower bearings on 
fully extended shock strut 
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1. = “ ial di * tan0<! between lower bearing and point o£ 

attachment of shock strut 
The seal friction force F is given by 

ffi 


s 

F = F . . .(2.71) 

" 1*1 

where 

u = coefficient of seal friction 

F = pneumatic force in the shock strut 

<a 

The total axial force in the shock strut is obtained as 


F = F + F + F, + F 

• a. h f j tm 


. . .(2.72) 


2.5. FORCES ON THE TYRE 

The exponential force-deflection of the tyre, 
neglecting hysteresis can be represented by 


F - m z* . . .(2.73) 

where 

F = vertical force applied to the tyre at the ground 

m , r = constants corresponding to various regimes of tyre 
deflection process 

z * radial deflection of tyre 

The drag load F, in the tyre ground interface is given 
by 

F » u F before wheel spin-up 

hg g 

s o 


after wheel spin-up 


. . .(2.74) 



where y is the coefficient of skid-friction between tyre and 
s 

ground 


The wheel spin-up time can be obtained from simple 
impulse momentum relationship. Thus , 


t 

au 


f y F dt 

J * v * 



. . . ( 2 . 75 ) 


where 



= moment of Inertia of the wheel assembly about 
its axle 

» horizontal velocity of landing of the aircraft 
* radius of the deflected shape of the tyre 
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CHAPTER 3 


RESULTS AND DISCUSSIONS 

The dynamic response of an aircraft during landing have 
been obtained. The equations of motion developed in chapter 2 
requires the numerical value of several quantities which may 
not be known to a designer before hand A parametric study is 
therefore carried out to determine the relative accuracy with 
which any parameter has to be known in advance in order to 
determine their Influence on the aircraft behaviour. 

The coefficient of orifice discharge plays the most 
significant role in any mathematical modeling of the landing 
gear. The dynamic loads on the aircraft can vary over 25% with 
the variations in C [6]. Due to lack of sufficient 

a 

experimental data on the air compression process in the shock 
strut, the exact value of the polytropic Index is also not 
known In advance. Hence these two parameters were varied to 
study their effect on the landing gear performance. The study 
also includes the effect of orifice diameter, initial air 
volume and initial air pressure. Besides the variation of the 
parameters of the shock strut for the heave model, the relative 


locations of 

the nose and main 

landing 

gears 

with 

respect 

to 

the aircraft 

center of gravity 

and the 

radius 

of 

1 

gyration 

of 


the aircraft for the heave-pitch model are also examined. 
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The equations of motion of the masses developed in 
sections 2. £2 and 2. <3 forms a set of coupled non-linear 
differential equations equal in number to the number of degrees 
of freedom of the airframe-landing gear system. These equations 
are solved by numerical integration. The method used for 
accomplishing this is the Gear method, vhich is a modification 
of Variable Order Variable Step Adam’s method. The equations of 
motion are reduced to a system of first order differential 
equations vhich are solved using a standard routine D02EBF in 
DEC 1090 mainframe. 

To start with, the programme reads in all the relevant 
physical characteristics of the aircraft, the landing gears and 
the initial conditions of the system. Using these informations 
the forces on the tyre and the geometry of the landing gear at 
any instance are determined. A distinction must be made at this 
point between the case when the force in the shock strut is 
greater than its pre-load and the case when it is less than its 
pre-load. In the latter case the force in the shock strut is 
treated as an unknown and an additional constraint is imposed 
on the system that there is no relative motion between the 
sprung and the unsprung masses as the shock strut is locked. 

The output information includes the displacement, 
velocity and acceleration of the system co-ordinates and the 
hydraulic, pneumatic and frictional forces In the shock strut. 
The output so obtained are plotted inorder to present the 
results in a graphical and tabular form for easy inference. 
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The aircraft and the landing gear data used 


for 


generating the results are placed in the Appendix. 


3.1 PARAMETRIC STUDY OF THE HEAVE MODEL 

3.1.1 Effect of Orifice Discharge Coefficient 

The influence of orifice discharge coefficient on the 

behaviour of the landing gear for the heave model are shown in 

figures 3.1a to 3.1k. Results are presented for a range of 

values of C from 0.2 to 1.0, when all other conditions are 
d 

assumed identical. The negative values of the displacement of 
the unsprung mass indicates an airborne state of the model. A 
decrease in C, increases the force m the shock strut due to a 

Ci 

higher hydraulic force (Table 1). Due to the larger force 
acting upwards on the sprung mass, the following observations 
can be justified. 

(a) . A decrease in the maximum downward displacement as 

decreases .( f igures 3.1a and 3.1b and Tables 2a and 2b) 

(b) .The aircraft ballooning back into air for the lower 

values of C, (figures 3.1c and 3. Id), the duration being 

d 

greater for C, = 0.2 compared to C, = 0.4. 

id d 

(c) .A higher maximum upward acceleration with a decrease in 

the value of C, (figures 3.1h and 3.1i and Table 2a and 

d 

2b) 

It can also be noted from the graphs that the downward 
acceleration of the sprung and unsprung masses are constant 
during the period the aircraft is fully airborne (figures 3.1h 
and 3 . 1 j ) . Therefore the variation of the velocity of the 
masses are linear during this time (figures 3 . le and 3. If). 
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PIG 3.1a VARIATION IN DISPLACEMENT OF THE SPRUNG MASS WITH 
ORIFICE DISCHARGE COEFFICIENT (0.2 & 0.4) 
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FIG 3.1b VARIATION IN DISPLACEMENT OF THE SPRUNG MASS WITH 
ORIFICE DISCHARGE COEFFICIENT (0.6, 0.8 & 1.0) 
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PIG 3.1c VARIATION IN DISPLACEMENT OF THE UNSPRUNG MASS WITH 
ORIFICE DISCHARGE COEFFICIENT (0.2 & 0.4) 
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FIG 3. Id VARIATION IN DISPLACEMENT OF THE UNSPRUNG MASS WITH 
ORIFICE DISCHARGE COEFFICIENT (0.6, 0.8 *1.0) 
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FIG 3. If VARIATION IN VELOCITY OF THE UNSPRUNG NASS WITH 
ORIFICE DISCHARGE COEFFICIENT (0.2 & 0.4) 
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PIG 3.1g VARIATION IN VELOCITY OF THE UNSPRUNG MASS WITH 
ORIFICE DISCHARGE COEFFICIENT (0.6, 0.8 & 1.0) 
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FIG 3.11 VARIATION IN ACCELERATION OF THE SPRUNG NASS WITH 
ORIFICE DISCHARGE COEFFICIENT (0.6, 0.8 & 1.0) 
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FIG 3.1k VARIATION IN ACCELERATION OF THE UNSPRUNG MASS WITH 
ORIFICE DISCHARGE COEFFICIENT (0. 6,0.8 & 1.0) 






Due to the higher force acting downwards on the unsprung 
mass, their peak downward displacements increase with 
decreasing (figures 3.1c and 3. Id). 

Table 1 gives the maximum values of pneumatic, hydraulic 
and the total force in the shock strut, the shock strut stroke 
and the force on the tyre. 


c d 

F 

a 

Newtons 

F k 

Newtons 

F 

m 

Newtons 

F 

t 

Newtons 

s 

cm 

0.2 

11345.0 

105110.0 

113600.0 

61492.0 

3.149 

0.4 

13514.0 

74056.0 

83679.0 

43992.0 

4.167 

0.6 

18836.0 

59290.0 

69989.0 

36123.0 

5.707 

0 . 8 

25817.0 

50112.0 

61633.0 

31459.0 

6.800 

1.0 

33438.0 

42436.0 

54548.0 

27534.0 

7.494 


TABLE 1. Variation in the peak values of the shock strut forces 
and tyre forces and the shock strut stroke with 

orifice discharge coefficient 

Because of the increased force in the shock strut with 

decreasing C , the ground reaction also increases with 
« 

deceasing orifice discharge coefficient. The rate of decrease 
of the tyre and the shock strut forces however is higher for 
lower values of C . Inorder to provide the higher pneumatic 

d 

force, the maximum strut closure increases with C, hence some 

d 

extra allowance should be made in the strut. for tfariations in 
Tables 2a and 2b shows the extreme positive and negative 
displacement and acceleration of the two masses for different 
values of C. . 

d 
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Maximum Displacement (cm} 

Maximum Acceleration (g) 

c d 

Up 

Down 


Up 

Down 

0.2 

15.5 

12.1 


7.2 

0.8 

0.4 

0.0 

13.4 


4.8 

0.8 

0.6 

• 

14.7 


3.9 

0.5 

O.S 

• 

16.2 


3.2 

0.4 

1.0 

* 

18.3 


2.7 

0.4 

TABLE 2a 

Variation in 

the extreme 

values 

of displacement and 


acceleration 

of the 

sprung 

ma.a& 

with orifice 


discharge coefficient. 





Maximum Displacement (cm) 

Maximum Acceleration (g) 

C d 

Up 

Down 

Up 

Down 

0.2 

16. S 

9.9 

13.2 

6.3 

0.4 

3.8 

8.9 

6.8 

3.4 

0.6 

. 

7.8 

8.6 

5.0 

0.8 

* 

6.8 

10.1 

4.0 

1.0 

- 

5.0 

11.5 

2.0 

TABLE 2b 

Variation in 

the extreme 

values of displacement and 


acceleration 

of the 

unsprung mass 

with orifice 


discharge coefficient. 
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Comparison o£ the figures and the values in Table 1 and 

Tables 2a and 2b show that magnitude of orifice coefficient 

has significant effect on the landing gear behaviour. The 

accurate prediction of the numerical value of C is therefore 

dl 

necessary for obtaining reliable results. 

3.1.2 Effect of Orifice Diameter 

A decrease in orifice diameter essentially Increases the 

hydraulic force in the shock strut (Table 3). Its effect is 

therefore similar to that of decreasing C, . Figures 3. 2. a to 

d 

3.2.£ and Tables 4a and 4b show the variation in the dynamic 
response of the aircraft for orifice diameters of 44 mm, 52 mm 
and 60mm. 

Table 3 gives the maximum values of pneumatic, hydraulic 
and the total force in the shock strut, the shock strut stroke 
and the force on the tyre. 


D 

n 

mm 

F 

OL 

Newtons 

F h 

Newtons 

F 

* 

Newtons 

F t 

Newtons 

s 

cm 

4.4 

14032.0 

71613.0 

81359.0 

42663.0 

4.365 

5.2 

18530.0 

59729.0 

70388.0 

36353. Q 

5.642 

6.0 

25265.0 

50741.0 

62194.0 

31742.0 

6.735 

TABLE 3 

Variation 

in the peak 

values of the shock strut forces 


and tyre 

forces and 

the shock 

strut stroke with 


orifice diameter 
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ORIFICE DIAMETER 



DISPLACEMENT Of UNSPRUNG MASS 



FIG 3.2b VARIATION IN DISPLACEMENT OF THE UNSPRUNG MASS WITH 
ORIFICE DIAMETER 
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FIG 3.2c VARIATION IN VELOCITY OF THE SPRUNG MASS WITH 
ORIFICE DIAMETER 
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FIG 3. 2d VARIATION IN VELOCITY OF THE UNSPRUNG MASS WITH 
ORIFICE DIAMETER 
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FIG 3.2© VARIATION IN ACCELERATION OF THE SPRUNG MASS WITH 
ORIFICE DIAMETER 
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Tables 4a and 4b shows the extreme positive and negative 
displacement and acceleration o£ the two masses for different 
values of D . 


Maximum 

Displacement (cm) 

Maximum 

Acceleration 

D n C»*0 

Up 

Down 

Up 

Down 

4.4 


13.5 

4.7 

0.8 

5.2 

- 

14.6 

3.9 

0.6 

6.0 


16.1 

3.2 

0.2 


TABLE 4a Variation in the extreme values of displacement and 
acceleration of the sprung mass with orifice 
diameter 


Naximum Displacement (cm) Maximum Acceleration (g) 
D^Cram) Up Down Up Down 


4.4 

3.1 

8.8 

7.0 

6.0 

5.2 

• 

7.8 

8.5 

5.4 

6.0 


6.9 

10.0 

4.1 


TABLE 4b Variation in the extreme values of displacement and 
acceleration of the unsprung mass with orifice 
diameter 
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3.1.3 Effect of Polytropic Index 

As the nature of the air-compression in the shock strut 
is not well defined, the value of n assumed varied from 1.0 
(isothermal) to 1 . 4 (adiabatic). Figures 3.3a to 3.3f show the 
influence of n on the aircraft response for the same set of 
initial conditions. The effect of increase in polytropic index 
is to increase the pneumatic force in the shock strut. For the 
present case, variation in air-compression process results in 
no significant difference in the velocity and acceleration 
response of the two masses. The sprung mass displacement shows 
some variation in the rebound after the initial impact, being 
lower for higher value of n (figure 3.3a). 


3.1.4 Effect of Initial Air-Volume 

V 

The parametric study for the initial air-volume is 

shown in figures 3.4a to 3.4f. The initial air-volume is fixed 

V 

at 275cc, 350cc and 425cc. A decrease in the value of a. 

o 

results in a higher pneumatic force for the same strut closure 
-acting upwards on the sprung mass. The displacement of the 
sprung mass shows some divergence in behaviour after the 
initial impact. Uhile the frequency and amplitude of 
oscillation do not have significant changes, their mean values 
are considerably different. This is due to the fact that oil is 
less compressible compared to air. Thus the static equilibrium 
position of the sprung mass is closer to the ground level a# 
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FIG 3.3b VARIATION IN DISPLACEMENT OF THE UNSPRUNG MASS WITH 
POLYTROPIC INDEX 
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FIG 3.3c VARIATION IN VELOCITY OF THE SPRUNG MASS WITH 
POLYTROPIC INDEX 
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FIG 3.3d VARIATION IN VELOCITY OF THE UNSPRUNG MASS UITH 
POLYTROPIC INDEX 
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FIG 3 . 3e VARIATION IN ACCELERATION OF THE SPRUNG NASS WITH 
POLYTROPIC INDEX 
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FIG 3.3£ VARIATION IN ACCELERATION OF THE UNSPRUNG MASS UITH 
POLYTROPIC INDEX 
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FIG 3.4a VARIATION IN DISPLACEMENT OF THE SPRUNG MASS WITH 
INITIAL AIR VOLUME 
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FIG 3.4b VARIATION IN DISPLACEMENT OF THE UNSPRUNG MASS WITH 
INITIAL AIR VOLUME 
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FIG 3 . 4e VARIATION IN ACCELERATION OF THE SPRUNG MASS WITH 
INITIAL AIR VOLUME 
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FIG 3 . 4f VARIATION IN ACCELERATION OF THE UNSPRUNG HASS WITH 
INITIAL AIR VOLUME 



increases (figure 3.4a) The peak acceleration of the 
unsprung mass also shows some variation. However on the whole 

it appears that for the present case the behaviour of the 
landing gear is not appreciably affected with variations in 
initial air volume. 


3.1.5 Effect of Initial Air - Pressure 

The influence of the initial air pressure on the 

response of the system are shown in figures 3.5a to 3.5f. 
Initial air pressures of 1 M Pa, 2 M Pa and 3 M Pa were used 
for generating system responses. An increased initial air 
pressure produces an Increase in the pneumatic force in the 
shock strut. Thus the sprung mass displaces more for lower 

p 

values of »* (figure 3.5a). As the aircraft is supported by 

the force in the shock strut at the static equilibrium position 

p 

the shock strut closure will be more for lower values of . 

Thus with decreasing initial air pressure, the sprung mass 
rests closer to the ground level (figure 3.5a). The velocity 
and acceleration responses of the sprung mass (figure 3.5c and 
3.5e) and the unsprung mass (3.5d and 3.5£) and the 
displacement of the unsprung mass (figure 3.5b) are not 

considerably affected by the initial air pressure In the 

present case, except that the amplitude of the responses is 
higher for lower ^ . In conclusion the effect of initial air 
pressure Is only marginal. 
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FIG 3.5b VARIATION IN DISPLACEMENT OF THE UNSPRUNG MASS WITH 
INITIAL AIR PRESSURE 
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FIG 3.5c VARIATION IN VELOCITY OF THE SPRUNG NASS WITH 
INITIAL AIR PRESSURE 
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FIG 3 . 5d VARIATION IN VELOCITY OF THE UNSPRUNG MASS WITH 
INITIAL AIR PRESSURE 
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FIG 3.5® VARIATION IN ACCELERATION OF THE SPRUNG MASS UITH 
INITIAL AIR PRESSURE 
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FIG 3 . 5f VARIATION IN ACCELERATION OF THE UNSPRUNG MASS WITH 
INITIAL AIR PRESSURE 



3.2 PARAMETRIC STUDY OF THE HEAVE-PITCH MODEL 


3.2.1 Effect of the Pitching Moment on Inertia of the Aircraft 
Figures 3.6a to 3.61 explains the effect of the pitching 
moment of inertia I of the aircraft for the same set of 
initial conditions . An increase in I y for the same mass of the 
aircraft produces a large inertia force in pitching. Thus, 

(a) . The displacement, velocity and acceleration of the 

airframe are not significantly influenced by I 

yy 

(figures 3.6a,3.6e and 3.6i) 

(b) . The pitching response of the aircraft is 

considerably affected by changes in the value of I 

yy 

(figures 3. 6b, 3. 6 f and 3.6j) 

(c) .Time history curves for the rigid body rotation show 

that the frequency of rotation Is less as I i s 

yy 

increased (figures 3.6b, 3.6f and 3.6 j) 

(d) . The response of the main wheels which are closer to 

the aircraft center of gravity are not significantly 
affected (figures 3.6d, 3.6h and 3.61) 

(e) .The nose wheel responses are appreciably affected as 

they are at some distance from the center of gravity 
of the aircraft (figures 3.6c, 3.6g and 3.6k) 

(f) .It can be concluded that the cockpit acceleration 

may also be considerably affected because of the 
long rotation arm 

Thus, while the effect of I yy may not be significant at 
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FIG 3 . 6d VARIATION IN DISPLACEMENT OF THE MAIN WHEEL WITH 
THE PITCHING MOMENT OF INERTIA 
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FIG 3 . 6e VARIATION IN VELOCITY OF THE SPRUNG MASS WITH 
THE PITCHING MOMENT OF INERTIA 
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FIG 3 . 6£ VARIATION IN VELOCITY OF PITCHING OF THE SPRUNG MASS 
WITH THE PITCHING MOMENT OF INERTIA 
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FIG 3.6g VARIATION IN VELOCITY OF THE NOSE WHEEL WITH 
THE PITCHING MOMENT OF INERTIA 
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FIG 3 . 6h VARIATION IN VELOCITY OF THE MAIN WHEEL WITH 
THE PITCHING MOMENT OF INERTIA 
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FIG 3 6} VARIATION IN ACCELEARTION OF PITCHING OF THE SPRUNG MASS 
WITH THE PITCHING MOMENT OF INERTIA 
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FIG 3.6k VARIATION IN ACCELERATION OF THE NOSE WHEEL WITH 
THE PITCHING MOMENT OF INERTIA 



FIG 3.61 VARIATION IN 
THE PITCHING 


locations close to the 


at points farther 
critical 


a ircraft center of gravity, its 
aircraft center of gravity 


ef f ect 
can be 


3-2*2 Ef f ect of banding t i . , . 

««ar Location with respect to the 

Aircraft center of Gravity 

The change* In aircraft dynMlc raaponaa ulth varlatll>na 
in landing gear layout is shown in figures 3.7a to 3.71. While 
the wheel-base is kept constant at 3.45 meters for the 
present study, three val Ues 40 Cfflf 56 Cffl and 7Q ^ £<>r the 

distance 1^ of the main wheel from the aircraft center of 
gravity were used in generating the response data. From the 
figures , 

(a) . Relative location of landing gears do not have 

considerable effect on the airframe response except 
at peak values. For the case studied, closer the 
location of the main gear to the aircraft center of 
gravity higher the peak displacement (3.7a), the 
peak velocity (3.7e), peak acceleration (3.71). The 
main wheel response is also found to be similar to 

that of the airframe (3.7d, 3.7h and 3.71) 

(b) . Displacement , velocity and acceleration response of 

the aircraft in pitching is observed to be higher 

for lower value of 1 (figures 3.7b, 3.7f and 3.7j) 

1*1 

(c) .For the lowest value of 1 , the nose wheel of the 

aircraft is found to leave the runway during rebound 
(figure 3.7c) 
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FIG 3.7b VARIATION IN PITCHING RESPONSE OF THE SPRUNG MASS WITH 
THE DISTANCE OF AIRCRAFT C.G FROM THE MAIN GEARS 
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DISPLACEMENT OP MAIN WHEEL 



FIG 3 . 7d VARIATION IN DISPLACEMENT OF THE MAIN WHEEL WITH 
THE DISTANCE OF AIRCRAFT C.G FROM THE MAIN GEARS 
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FIG 3 . 7e VARIATION IN VELOCITY OF THE SPRUNG MASS WITH 

THE DISTANCE OF AIRCRAFT C.G FROM THE MAIN GEARS 
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FIG 3 . 7£ VARIATION IN VELOCITY OF PITCHING OF THE SPRUNG NASS 
WITH THE DISTANCE OF AIRCRAFT C.G FRON THE MAIN GEARS 



VELOCITY Or MAIH TO EEL 



< S/TW 5 'P’S 


1 04 - 


FIG 3 . 7h VARIATION IN VELOCITY OF THE MAIN WHEEL WITH 

THE DISTANCE OF AIRCRAFT C.G FROM THE MAIN GEARS 





FIG 3 , 7i VARIATION IN ACCELERATION OF THE SPRUNG HASS WITH 
THE DISTANCE OF AIRCRAFT C.G FROH THE HAIN GEARS 




FIG 3 . 7 j VARIATION IN ACCELEARTION OF PITCHING OF THE SPRUNG MASS 
WITH THE DISTANCE OF AIRCRAFT C.G FROM THE MAIN GEARS 
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(d) -Considerable difference in the peak 

acceleration of the main and nose wheel can be 
observed with changes in 1 (figures 3.7k and 3.71) 
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Tables 5a, 5b and 5c 

gives 

the peak 

values 

of 

displacement and acceleration of 

the 

sprung mass, 

the 

nose 

wheel and main wheel with variations in 

the distance 

of 

the 


main gear from the aircraft center of gravity. 


Maximum Displacements 

! a/c c.g pitching nose-wheel main-wheel 

m 


cm 

cm 

rad 

cm 

cm 

40.0 

30.9 

0.118 

5.4 

8.6 

56.0 

30.0 

0.074 

4.8 

8.2 

70.0 

28.8 

0.053 

4.8 

7.5 

TABLE 5a 

Variation in 

the distance 

the extreme values of 

of the main gear 

displacement with 

from the aircraft 


center of gravity 
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Maximum Acceleration of the airframe 

<-» 

vertical 

translation 

pitching 

00 

6 

vw 

O 


positive 

negative 

positive 

negative 

cm 

8 

8 

rad/s .a 

rad/s. s 

40.0 

0.68 

2.4 

37.0 

23.0 

56.0 

0.60 

2.3 

23.0 

18.0 

70.0 

0.58 

2.1 

19.0 

14.0 

TABLE 5b 

Variation in 

the extreme 

values of acceleration 


o£ the airframe with the distance of the main gear 
from the aircraft center of gravity 


Maximum Acceleration of the wheels 


L nt 

Nose 

wheel 

Main 

wheel 


positive 

negative 

positive 

negative 

cm 

8 

8 

8 

8 

40.0 

11.0 

14.0 

2.2 

19.5 

56.0 

4.9 

11.0 

2.1 

8.8 

70.0 

2.9 

10.0 

l. i 

8.4 


TABLE 5c Variation in the extreme values of acceleration 

of the wheels with the distance of the main gear from 
the aircraft center of gravity 
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Thus even for the same mass of the aircraft and for the 
same initial conditions the relative location of the main and 
nose gears with respect to the aircraft center of gravity has a 
considerable influence on the dynamic response of the aircraft. 
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CHAPTER Ar 


CONCLUSIONS AND RECOMMENDATIONS 

An analysis of a two degrees of freedom heave model of 
an articulated landing gear has been presented The variation 
in the dynamic response of the system to changes in the orifice 
discharge coefficient, polytropic index of air compression 
process, diameter of the orifice, initial air volume and 
initial air pressure have been studied. 

The study on the heave model has been extended to a 
heave-pitch model of an aircraft with articulated nose landing 
gear and telescopic main gears. The effect of the radius of 
gyration and the landing gear layout configuration also has 
been included. 

4.1 CONCLUSIONS 

The main conclusions obtained regarding the landing gear 
behaviour at touch-down from this study on the heave and the 
heave-pitch models are presented below. 

4.1.1 The Heave model 

The conclusions drawn from the study of this model are : 

1. The aircraft response is strongly affected by the orifice 
diameter and the coefficient of discharge. The changes in 
their values influence the response characteristics in a 
similar fashion. The peak displacement response increases 
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with increase in the values of these parameters. 

2. Even though the thermodynamic law of air compression 
process in the shock strut is not well known, it has only 
a marginal influence on the dynamic response of the 
aircraft . 

3. Initial air volume and air pressure do not have 
appreciable effect on the aircraft response. However the 
static equilibrium position of the airframe is affected to 
some extent. 

4. The performance parameters for a landing gear are usually 
difficult to evaluate accurately. Some extra allowance 
should be made in the shock strut to take into account the 
variable nature of some of these parameters. 

4.1.2 The Heave - Pitch Model 

The study of the heave-pitch model helps in drawing the 
following conclusions : 

1. Even for the same mass of the aircraft its relative 

distribution plays an important part in landing 

performance especially at locations farther from the 
aircraft center of gravity. 

2. The layout of the landing gear with the aircraft center of 
gravity significantly alters the peak values of response. 
Landing gear locations may be optimized for improving 
performance. 
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4.2 SUGGESTIONS FOR FURTHER INVESTIGATION 

The following suggestions are made for further study in 
this area. 

1. A more realistic estimate of aircraft response can be 
obtained by accounting for the flexibility of airframe. 

2. In practice the landing may not be perfectly symmetric. 
The effect of aircraft roll and that of the cross wind on 
the airframe may be considerable and should be studied to 
obtain the maximum landing stresses. 

3. The Influence of runway unevenness on the aircraft during 
its ground roll in landing can be studied and their effect 
on structural fatigue and fatigue of undercarriage load 
attachment points determined. 

4. Pilot’s input to the system at touch-down can also be 
included. 
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APPEND ix 


AIRCRAFT DATA 


i . general 

1 l 


1. Gross weight of the aircraft 

2. Height of the stain wheels 

3. Height of the nose wheels 

4. Uheel base 

5. Distance of aircraft c.g from 
main wheels 

6. Distance of aircraft c.g from 
nose wheel 

7. Distance of aerodynamic center 
from aircraft c.g 

8. Moment of inertia of the 
aircraft in pitching 

9. Landing velocity of aircraft 


= 40750.7 Newtons 
= 627.0 Newtons 
= 55.0 Newtons 
= 3.45 m 
= 0.56 m 

= 2.89 m 

= 0.2m 

= 1000 kg-m 2 
= 54 m/s 


1 1 . CHARECTERI STI CS OF THE SHOCK STRUT OF THE MAIM GEAR 


1. Coefficient of discharge 

\ i 

2. Coef f icient of seal friction 

3. Mass density of field 
4. Orifice area _ j' 


if 
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* 0.6 
* 0 - 

= 850 kg/cu.m 
= 2.405 sq.cm 


5 . Hydraulic area 

6 . Pneumatic area 

7 . Extended air pressure 

8. Ext ended air volume 

9. Index for air compression 
10. Initial oil volume 

11. Full stroke 

12. Coefficient of friction at 
bearings 

13. Axial distance between bearings 

14. Coulomb friction 

15. Preload 

III . CHARECTERI STICS OF THE SHOCK STRUT 


1. Coef f icient of discharge 

2 . Coefficient of seal friction 

3. Hass density of fluid 
4. Orifice area 

5. Hydraulic area 

6. Pneumatic area 

7. Extended air pressure 

8. Extended air volume 

9. Index for air compression 
18. Initial oil volume 
11. Full stroks 
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s 63.6 sq . cm 
* 49.5 sq.cm 
■ 140 psi 
= 1462.0 c.c 
= 1.3 

= 1950.0 c.c 
= 26.4 cm 
= 0.05 to 0.1 

= 20.0 cm 
= Q 

= 4780.8 Newtons 

OF THE NOSE GEAR 


= 0.6 
* 8 

* 850 kg/cu.m 
= 0.216 sq.cm 

* 42.0 sq.cm 
= 27.2 sq.cm 

* 400 psi 

= 273.6 c.c 
= 1.3 

= 400.0 c.c 
■ 9.3 cm 



12. Coefficient of friction at 
bearings 

13. Axial distance between bearings 

14. Coulomb friction 

15 . Preload 

IV. MAIN TYRES 


1. Radius 
2. Spin-up time 

3. Polar moment of inertia 

4. Load-deflection curve given 

V.NOSE TYRE 


1 .Radius 
2. Spin-up time 

3. Polar aonent of inertia 

4 . Load-deflection curve given 

V. DIMENSIONS OF THE MOSE GEAR ( refer 


Aj s j.O c« 

A 2 * 20.0 ca 

A 3 * 40.0 a 

A. = 3.0 cn 

4 


= 0.05 to 0.1 

= 5.0 cm 
- G 

= 7512.4 Newtons 


a 25.0 cm 
= 0.038 secs 
= 3217 kg-cm 2 


= 20.0 cm 
= 0.030 secs 
=1173 kg-cm 2 

fig 2.1 ) 


. Ilf 



50.0 ca 
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